Dynamical description of vesicle growth and shape change 
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We systematize and extend the description of vesicle growth and shape change using linear 
nonequilibrium thermodynamics. By restricting the study to shape changes from spheres to ax- 
isymmetric ellipsoids, we are able to give a consistent formulation which includes the lateral tension 
of the vesicle membrane. This allows us to generalize and correct a previous calculation. Our present 
calculations suggest that, for small growing vesicles, a prolate ellipsoidal shape should be favored 
over oblate ellipsoids, whereas for large growing vesicles oblates should be favored over prolates. The 
validity of this prediction is examined in the light of the various assumptions made in its derivation. 
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INTRODUCTION 



Vesicles are fascinating structures for several rea- 
sons [l|: they occur in a wide range of shapes and sizes, 
they are cell-like and are frequently used in the modeling 
of proto-cells and they present a formidable task for the 
mathematical modeler. It is not hard to see why they 
are so difficult to describe theoretically. They are in ef- 
fect a closed membrane made up of a lipid bilayer, which 
allows water and solutes to permeate from the environ- 
ment into the interior of the vesicle While they are 
frequently modeled as a two-dimensional surface, the fi- 
nite thickness of the bilayer as well as its structure, plays 
an important role in understanding their behavior [2| . In 
a pioneering work, Helfrich 0, H| used the analogy be- 
tween the rod-like shape of nematic liquid crystals and 
the lipids forming the vesicle bilayer, to write down an 
expression for the energy of the membrane. He envisaged 
the membrane as a curved surface, and accounted for the 
fact that the membrane had structure through the intro- 
duction of a phenomenological constant, Co, the spon- 
taneous curvature. Since then a variety of other models 
have been proposed, most of which start from the idea of 
a purely geometrical surface to represent the membrane, 
with additional structure introduced in various ways 

For most of the period since these models were pro- 
posed the main focus of their study has been to look for 
the shapes with the smallest energy under fixed condi- 
tions of constant volume, V and constant surface area, A. 
It was argued that this would give the expected shape of 
the vesicle at these values of V and A, and indeed a range 
of shapes emerged from this analysis 0] ■ However this is 
a purely static approach — there is no mechanism which 
postulates how the transitions between different shapes 
occurred, or how long these changes take. We recently 
carried out a preliminary study of a dynamical process 
designed to describe the change in vesicle shape Q ■ The 
change in surface area occurred because of the slow accre- 
tion of lipids onto the surface from the fluid surrounding 
the vesicle. This resulted in a change in the volume of 
the vesicle due to the influx of water and solutes through 



the membrane. The process was assumed to be suffi- 
ciently slow that the formalism of linear nonequilibrium 
thermodynamics (LNET) could be used [6[. The expec- 
tation was to find stability conditions which determine 
the shape of the vesicle at various stages of its growth. 

While we believe that modeling vesicles as surfaces 
with an energy of the Helfrich type, which change their 
shape and size according to LNET, is a minimalist de- 
scription which is capable of answering many interesting 
questions, it is still extremely mathematically complex. 
The surface may be of an arbitrary shape, and so requires 
the formalism of differential geometry for its study, the 
actual dynamics, as opposed to the statics, of mixtures 
of the rod-like lipids and the point-like fluid molecules is 
very non-trivial, and the nonequilibrium thermodynam- 
ics of discontinuous structures such as vesicles as been 
little studied because of its difficulty. Some of these ques- 
tions will be discussed in a forthcoming publication Q, 
where we will develop the formalism in a systematic fash- 
ion. However in our view it is also valuable to proceed in 
stages, and build on what has gone before [5lla-ITl|. and 
not to attempt to introduce all aspects of the problem si- 
multaneously. This is the philosophy that we will follow 
in this paper; we will extend our previous calculation to 
include surface tension (frequently called lateral tension 
when discussing vesicles) in the context of axisymmetric 
ellipsoidal shapes, treating the area A, as well as the vol- 
ume, V , as thermodynamic variables. This will form the 
basis for the treatment of more general shapes in the fu- 
ture. We will also correct some of the analysis presented 
in @. 



II. FORMALISM 

We shall in common with most other authors adopt 
the spontaneous curvature model of Helfrich Q in which 
the energy of the membrane is given by 



E m = | £ (2H - Co) 2 dA. 



(1) 
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where H is the local mean curvature, dA is an element of 
the surface A and k is the bending rigidity. Since a purely 
static analysis of Eq. (JT]) shows that when a spherical sur- 
face becomes unstable it is replaced with an ellipsoid [l2j , 
we will restrict ourselves to shapes which are spheres or 
axisymmetric ellipsoids. These will be parametrized in 
Cartesian coordinates by 

x = a sin 9 cos 4> , 
y = a sin 9 sin <f> , 

z = c cos 6 , (2) 

where < (f> < 27r, < 9 < -k and where a and c are 
constants. For a sphere a = c = r, the radius. We 
have not included a term proportional to the Gaussian 
curvature in the expression for the energy, since this is a 
constant for the kinds of surfaces that we are considering 
here 

The surface area, volume and membrane energy ((T|) can 
all be found in closed form for an axisymmetric ellipsoid 
and are given explicitly in the appendix of [B| (there is 
however a typographical error in Eqs. (A9) of that paper: 
the minus sign in the factor — 2a 3 /3c 4 should be absent). 
They are all functions of the two variables a and c which 
characterize the ellipsoid, with E m additionally depend- 
ing on the parameters k and Cq. Since the volume takes 
the simple form V = 47ra 2 c/3, c may be explicitly elim- 
inated in favor of a and V in the expressions for A and 
E m . Subsequently A = A(V,a) may be inverted to give 
a = a(V, A). In this way we see that the energy of the 
membrane is a function of V and A only. Since in the 
Helfrich model the membrane is completely described by 
this energy, we deduce that the thermodynamic descrip- 
tion of the membrane can be achieved using these two 
variables. 

As previously indicated we use LNET to describe the 
time evolution of the vesicle. Throughout, we assume 
that the temperature, T, is fixed, since we are not inter- 
ested in changes in shape which come about because of 
a change in temperature. For simplicity, we assume that 
no solutes are present, although it is not too difficult to 
extend the treatment to include them. 

The system is assumed to comprise two regions, the in- 
terior and exterior, that are separated by a third region, 
the membrane. This third region is assumed to be very 
thin, and our aim is to modify the thermodynamic de- 
scription so that the membrane can (following Helfrich) 
be simply regarded as a boundary, or geometrical surface, 
between the first two regions. In the formalism of LNET, 
each region is assumed to be in local equilibrium @ , with 
the thermodynamic relation TdS = dE — dW holding in 
each region. Here S is the entropy, E the internal en- 
ergy and W the work done on the system [l3| . The three 
regions will be labeled i, e and m respectively. 

The internal energy and work done have the following 
forms: 

(i) The internal energy of the interior, Ei, and exterior, 
E e , will be that of the fluid in these two regions, 



in our case water. Their sum is denoted by E w . 
The internal energy of the membrane is given by 
Eq. CP). 

(ii) Considering the interior and exterior regions in- 
dividually, each will have work done on them of 
— PdV if their volumes change, where P is the 
pressure in that region. Adding these gives the 
total contribution for work of this type to be 
— PidVi — P e dV e . The membrane will be assumed 
to have negligible volume, and so gives no contri- 
bution. This also means that dVi = —dV e = dV, 
and so the total work done on the system due to the 
interior increasing its volume by dV is (P e — Pi)dV . 

(iii) There will also be a work done if the area of the 
membrane increases by dA, equal to adA, where a 
is the surface tension This would exist even 
if the system consisted of two regions of different 
fluids, with no membrane separating them. The 
existence of a membrane with non-trivial structure 
separating the interior and exterior, means that a 
will have a more complex functional form which 
reflects this structure. For this reason we will follow 
the usual usage in this field and refer to it as the 
lateral tension. 

Adding up the contributions from the three regions one 
finds 

TdS* = dE w + dE m - (AP) dV - a dA, (3) 

where now S is the total entropy of the system and AP = 
(P e — Pi) is the pressure difference between the exterior 
and interior. The membrane does not explicitly appear 
in the terms relating to work. We can also eliminate it 
from the internal energy by noting that, since E m is a 
function of V and A, 

This allows us to write the thermodynamic relation ([3]) 
for the system as 

TdS = dE w - (AP) cff d^ - ((7 eff ) dA, (5) 



where 




Equation ([5]) is the thermodynamic relation for two re- 
gions separated by a boundary with no material proper- 
ties. The effect of the membrane simply changes the pres- 
sure difference and lateral tension from AP to (AP) c g 
and from a to cr G ff respectively. Therefore, as long as we 
make these replacements, we may ignore the membrane 
from a thermodynamic point of view, and simply treat 
it as a boundary which separates the inside of the vesi- 
cle from the environment. In our previous analysis § we 
did not treat the area as an independent variable. This is 
clearly consistent for spherical vesicles, but not for those 
which have an axisymmetric ellipsoidal shape. 
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III. DYNAMICS 

One of the central features of LNET is the relation 
between the forces, X a , which cause the state of the sys- 
tem to change, and the fluxes, J a , which are the result of 
these changes @ . Within the formalism of LNET these 
are linearly related: J a = J^b LabXb, where the L a b are 
constants, the Onsager coefficients. The forces and fluxes 
can be identified in a systematic way Q , but for the spe- 
cific problem of interest to us here, where the relevant 
thermodynamic variables are simply V and A, we may 
proceed more directly. We will also restrict ourselves to 
just one of the fluxes: that due to water flowing through 
the membrane into the interior of the vesicle. This will be 
denoted by J w . The direct effect that causes this flux of 
water is the pressure difference between the exterior and 
interior regions [7JI, [lj| , which incorporating the effect of 
the membrane is (AP) g. 

If we only took into account this direct effect, as 
we did in our previous treatment [5|, we would write 
J w = Lp(AP) ff', where L p is the hydraulic conductivity 
of the membrane. However there will also be an indirect 
effect 0] for which the driver will be <7 c ff, and so 



Jw = Lp(AP) c ff + i CT CT off , 



(7) 



where L a is a second Onsager coefficient. We shall jus- 
tify the choices of forces and fluxes in more detail else- 
where , but we can give a simple microscopic argument 
showing how the flux proportional to the lateral tension 
comes about. For a positively curved membrane, the 
lipids are arranged in a funnel-like configuration which 
inhibits the flow of water across the membrane. Decreas- 
ing the curvature of a membrane at any point therefore 
permits a greater flow of water due to the lipids becom- 
ing increasingly parallel. In this way, for fixed shapes, an 
increase in the surface area causes a slight alignment of 
adjacent lipids and hence permits a greater flow across 
the membrane. This resultant flow is proportional to the 
lateral tension, the measure of how the membrane energy 
changes with area. 

We are now in a position to introduce the dynamics 
explicitly. The rate of increase of the volume of the vesicle 
will be given by dV/dt — AJ W , and so from Eqs. © and 
([7]) we have 



dt 



L, 



AP 



(dE n 



V dV 



fdE n 



V OA 



(8) 

We also need to describe how the surface area of the 
vesicle grows due to the inclusion of lipids from the envi- 
ronment. The simplest assumption is that these attach 
themselves uniformly over the entire surface, so that the 
surface grows at a constant rate per unit area which we 
will denote by A [1Q : 

d A 

— = XA => A(t) = A(0)e xt . (9) 

Before we go on to investigate the dynamics for an 
ellipsoidal shape, let us briefly consider the result for a 



sphere. In this case V and A are not independent vari- 
ables, and so the thermodynamic relation ([5]) should read 



TdS = dE w - (AP) s c p s hclc dV, 



(10) 



where 



(Ap) sphere ^ = {Ap) ^ ^ + ^ ^ 



APdV + adA 



dE„ 



dA 



= APdV + adA- dE„ 



AP 



f dE„ 



I dV 



dA 
dV 



dA 



dV, 



fdE n 



V dV 



dV 



(11) 



where in the last line the derivatives are not partial 
derivatives since for a sphere E m = E m (V) and A = 
A(V). We may calculate these in a straightforward way. 
The mean curvature of a sphere of radius r is and 
so from Eq. (pj 



E m = 2ttk (C r - 2Y 



dgm C K 



(12) 

Using dA/dV = 2/r, we find from Eqs. dTTJ) and (fT2"|) 

that 



(AP) 



sphere 



= AP — 



C()K 



2(7 



* --p-(C r-2) + T 



(13) 



However for a sphere J w — L p (AP) c s, and so in equilib 
rium when there is no flow of water, J w = 0, Eq. |T~ 
gives 



C K 



2a 



AP \- (Cor — 2) H = 



(14) 



This is a standard result from the studies of spherical 
vesicles in equilibrium [TJJ , if we make allowance for the 
different sign conventions for the pressure difference and 
the lateral tension. In previous work [jj [ll|, we did 
not include the lateral tension in our description, and 
so our equilibrium result did not include the final term 
in Eq. (IT4l . 

It should be emphasized that Eq. (|T4"]) is a consequence 
of asking that the vesicle is in static equilibrium, so that 
in particular no lipids are being added to the exterior 
surface, leading to no increase in the surface area A. A 
condition for dynamic equilibrium can also be obtained. 
This is a stationary state in which the vesicle remains 
turgid and grows like a sphere. In this case, recalling that 
the area grows according to Eqs. ©, one gets dV/dt = 
27rAr 3 . The spherical version of Eqs. © and ©, dV/dt = 
AL p (AP)^ horo , therefore gives 



Xr 
2L V 



AP 



— 2~ W° r " 



2a 



(15) 



a result that could not be derived from a purely static 
description. Equation (fT5"|) is Eq. (21) of [Bj], but with the 
lateral tension now taken into account. 
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The inclusion of the lateral tension is even more im- 
portant for axisymmetric shapes of the kind we are con- 
sidering here, because now A is an independent variable. 
Since we are concerned with questions of stability, we 
will only consider ellipsoids which differ in shape from 
the sphere very slightly. In this case the parameters a 
and c which describe the ellipsoid (see Eq. may be 
expressed as 



R(l + one) , c = R(l + cie) , 



(16) 



where e is a small quantity and a\ and c\ are numbers 
which characterize the shape of the ellipsoid: if a% > c\ 
it is oblate and if a\ < c\ it is prolate. The quantity R 
reduces to the radius of the sphere as e — > 0, but care is 
required in its definition, as we will discuss in more detail 
below. 

Using standard results and Eq. ([To]) , it is straight- 
forward to calculate the surface area and volume of the 
ellipsoid for small e. From the explicit forms given in the 
appendix of 5] it is found that 



.4 



AttR 2 



l + -(2a 1 +c 1 )e + 0(e 2 ) 



V = -ttR 3 [1 + (2oi + ci) e + O (e 2 )] . (17) 

However, as we will now show, it is not consistent to 
assume that R is independent of e if we assume a growth 
law of the form (|9]). To see this, we write the expression 
for A given in Eqs. (fl~7| as A — Air R 2 cf>(e) , where </>(e) is 
the expression in the square brackets. Then, 



XA = — = 
dt 

=> A = 



dt 



which implies that R(t) = e At / 2 [(^(i)] 1 ^ 2 , up to an over- 
all multiplicative constant. Since <f> is a function of e, so is 
R. Therefore, for consistency, we cannot use 7? when car- 
rying out a perturbative expansion in e, since it contains 
hidden e dependence. Instead we should use the radial 

1/2 

variable r(t) — R(t) [4>(t)] ' which is e— independent. 
This is equivalent to determining r through the condi- 
tion A = Airr 2 , for any axisymmetric ellipsoid. Clearly 
r is the radius of the sphere which has the same surface 
area as the ellipsoid. 

The correct procedure to investigate the dynamics 
of vesicle growth perturbatively in e is, therefore, to 
use closed form expressions from [5[ and Eq. (|16p . to 
determine the results (|17[) to the required order, but 
then to set A = 4irr 2 . This can be inverted to find 

]/2 

R(t) = r(t) [(f>(t)] 1 , allowing V and E rn to be found 
as functions of r and e. To first order, this procedure 
gives 



(18) 



1 + 3 (2a x 



ci)+0(e 2 



(19) 



leading to V = (47rr 3 /3) [l + C(e 2 )] . In fact, since A = 
Airr 2 exactly, the volume is directly related to the so- 
called reduced volume, defined by v — 6^/nV/A 3 / 2 , by 



47rr 3 
V = — - — v. 



(20) 



Although V has no term of order e when expressed in 
terms of the variable r, it does turn out to have a term 
of order e 2 : 



1 
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n2 2 
Ci) € 



0(e 3 ) 



(21) 



As a check we note that v < 1 for all cases except the 
sphere (ai = ci), as required. In our previous calcula- 
tion 0] we used the variable R, rather than r, which gave 
incorrect results for the coefficients in the expansion in 
e. Correcting these by eliminating R in favor of r, makes 
previously cumbersome results look far more elegant. For 
example, the energy of the membrane becomes 

E m (r, e) = £(°>(r) + a 2 E^(r)e 2 + O (e 3 ) , (22) 



where 



8ttk 



(r) = 2ttk (C r - 2) 2 , E^ (r) = — (6 - C r) , 

(23) 

and where 



4 2 
a 2 = TT («i - ci) • 
15 



(24) 



A purely static analysis of the stability of an ellipsoidal 
vesicle would compare the energy of the membrane (|22|) 
to the energy of a sp herical vesicle — the same equation, 
but with e = [12J. The conclusion would be that the 
ellipsoid is the more stable if E^(r) < 0, i.e. if C r > 6. 
However there is no dynamics in this picture at all. To 
achieve a more physical description of the time evolu- 
tion of the vesicle we utilize Eq. (|SJ). To do this we first 
need to evaluate (dE m / '8V)a and (dE m /dA)v, but we 
cannot proceed directly, since we know E m — E m (r, e) 
rather than E m = E m (V, A). To circumvent this prob- 
lem, suppose that we have inverted (in principle, not in 
practice) V — V(r,e) to obtain e = e(r,V). Then since 
E m = E m (r,e(r,V)), 



dE m 
dV 

dE m 
OA 

1 

87rr 



dE n 



(dE n 



dV ) r V 9e 

= J_ fdErn\ 

v 8nr V dr )v 
dE m \ ( dE m \ ( de 
dr 



8e_ 



dr 



(25) 



So we need only to calculate (de/dV) r and (de/dr)v- 
These may be found from Eq. (|2"U1) by noting that the 
reduced volume is a function only of e. Then 



V '^[dVj r 4vrr 3 /3' 



V 



r 47rr 3 /3' 
(26) 
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where v'(e) = dv/de. Substituting expressions (|26l) into 
Eqs. (1211) gives 



(dE n 



K(e)]" 1 (dE, 



V dV /A 
( dE m 
\ dA 



47rr 3 /3 \ de 



1 


fdE m " 


\ 3V 


(dE m \ 


87IT 


V dr , 


)~ 2A 


V dV J 



.(27) 



These two equations allow us to find (dE m /dV)A and 
(dE m /dA)v if we know E m and V as functions of r and 
e. 

Using the results ([20 |) -([24 |) to second order in e, the 
partial derivatives (|27|) are given by 



dE m 
dV 

dE m 
dA 



1 



1 

87rr 



47rr 3 /3 



£ (2) (r)+0(e) 



dr 



-E( 2 )(r 



■0(e) (28) 



It is important to note that while we have used the re- 
sults for E m and V correct to second order in e, we have 
only been able to calculate the required derivatives to ze- 
roth order in e. The reason for this can be traced back 
to v'(e) being of order e and the differentiation of E m 
with respect to e also giving an expression of order e. 
Together these reduce the powers of e by 2 in the calcu- 
lation of (dE m /dV)A- The result can however be used 
to check that we recover the previously derived form for 
the spherical vesicle in the limit e — > 0. Substituting 
Eqs. (EHD into Eq. ©, and using dV/dt = 27rAr 3 , valid 
for a sphere, one finds 



y =L P <AP 



5 (6 " Cor) 



This is identical to Eq. (|15p . already derived for the 
sphere provided that we make the identification L a = 
2L p /r. It should be emphasized that this identification 
is only being made to obtain the correct result in the 
spherical limit, bearing in mind that L a is not defined in 
this case. 



IV. STABILITY 

To determine if and when the shape of the vesicle 
starts to deviate from that of a sphere, we calculate the 
time derivative of the reduced volume. This will be zero 
(v = 1) when the vesicle remains spherical, but will start 
to decrease as soon as it adopts another shape. From 
Eq. (HOD, 



dV , , dr Airr 3 dv 

—— = 4irr v — — I — . 

dt dt 3 dt 



(30) 



The zeroth-order part of the first term on the right-hand 
side of Eq. (|3T)]) has already been included in the station- 
ary condition (|29[) . At first order in e, the second term 



in this equation is the relevant one. To find the right- 
hand side of Eq. flSJ) to order e it is necessary to find 
(dE m / '8V)a to order e 3 . This can be carried out in the 
same way as described above for the calculation to order 
e 2 , although the intermediate steps are sufficiently com- 
plicated that we used Mathematica [13] • Nevertheless, 
the final result is quite simple: 



x (5C r + 6) [2L P - rL< 



or defining a critical radius, r. 
de 



by r c 



.] e + 0(e 2 )(31) 
2L p / L a , 



a2 dt 



2nL a 
7r 4 



(ai -ci)(5CV + 6)(r-r c ) + 0(e) 



(32) 

This stability condition differs from the one we had 
derived in our earlier study of this question [||, in that 
it predicts a linear growth away from the sphere, and 
not an exponential one. Since a.2 > 0, if the right-hand 
side of Eq. (|32|) is positive the sphere is linearly unstable, 
if it is negative it is linearly stable. We see that oblate 
(ai > Ci) perturbations will destabilize a sphere if r > r c , 
whereas prolate (ai < c\) perturbations will destabilize 
a sphere if r < r c . The picture we have is that as the 
sphere grows, it is susceptible to fluctuations which give 
it a prolate ellipsoidal shape if it has a radius less than 
r c , and susceptible to fluctuations which give it an oblate 
ellipsoidal shape if it has a radius greater than r c . 

How do these predictions compare with the currently 
available experimental results? According to our predic- 
tions, spherical shapes should be unstable in a purely dy- 
namical setting where the membrane surface grows due to 
the successive inclusion of lipid molecules and where the 
volume then increases due to the inflow of water. Under 
static conditions, however, the Helfrich energy functional 
implies that spherical vesicles are stable below a critical 
radius. It can be argued that the reason why spherical 
shapes are often seen in experiments [lj is a reflection of 
the specific experimental setting adopted. If the dynami- 
cal mechanisms of the type considered here are essentially 
negligible, then static effects are expected to prevail and, 
consequently, spheres are the energetically favored con- 
figurations. 

In [l8j], the process of vesicle formation was investi- 
gated in a rich and dynamic phospholipid mixture. Dy- 
namical light scattering and transmission electron mi- 
croscopy experiments were performed to resolve the vesi- 
cles' shape and quantify their associated sizes. The ob- 
served vesicles were relatively large (hydrodynamic radii 
~ 200 nm or > 500 nm) and corresponded to either 
oblate ellipsoids or triaxial ellipsoids. The fact that a 
significant fraction of the (giant) vesicles belonged to the 
oblate ellipsoids family is consistent with the prediction 
from our analysis. 

Clearly, there are number of other effects (e.g. tem- 
perature fluctuations) which are implicated in the above 
mentioned experiments and not included in our modeling 
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efforts. We also need to consider more generic shapes, 
including triaxial ellipsoids. The probability that per- 
turbations are axisymmetric are presumably small, nev- 
ertheless, oblate ellipsoids are observed in experiments 
where the dynamics plays a role and the vesicles are suf- 
ficiently large. This is for instance the case in [l9|, [2(| , 
where initial spherical vesicles are experimentally shown 
to deform into oblate ellipsoids, rather then prolate ones, 
in the presence of osmotic driving pressures. This is a 
distinctive feature which is not captured by any static 
approaches to the problem of the morphology of vesicles 
and one that could perhaps be explained by a crucial 
interplay between shape and dynamics. 

V. CONCLUSIONS 

In this paper we have extended and corrected earlier 
work on the growth and stability of vesicles. The errors 
in our previous paper [5| were firstly an inappropriate 
choice of characteristic size for the ellipsoidal vesicle (i? 
rather than r). This was a subtle point relating to the 
fact that assuming the radius to be independent of the de- 
gree of deformation e was not compatible with the growth 
law we chose. The second error was one of omission: we 
did not include the lateral tension in our original anal- 
ysis. Finally, there was a simple typographical error in 
Eqs. (A9) of 0. While correcting the first point would 
not have changed the conclusions to any great extent, 
adding the surface tension does: the prediction is now 
that smaller vesicles should tend to be prolate and larger 
vesicles should tend to be oblate. 

There are several caveats that we should make in re- 
gard to this prediction. Foremost among these is that the 
model we have adopted is very simple, and even though 



it is the one frequently used in the literature, we should 
keep this constantly in mind. It should be possible to 
extend these results to the ADE [2l[ model, which is 
slightly more realistic. It would be interesting to see if 
the predictions are changed in any way. Another caveat 
is that we have assumed that the sphere becomes unsta- 
ble to an ellipse defined by Eqs. ©. Although this is 
what is found from investigations using variational tech- 
niques [12j . it would be more consistent not to assign a 
specific functional form to the shape of the membrane. 
This is currently being investigated; the analysis is con- 
siderably more complicated and the results will be pre- 
sented elsewhere [Tj. Finally, LNET for this problem is 
very underdeveloped, and rather complex, with plenty 
of scope for pitfalls. This will also be more extensively 
discussed elsewhere [7[. 

Even given all these caveats, we still believe that the 
current work is a significant step forward in understand- 
ing vesicle growth and shape changes. Virtually all pre- 
vious work was static and made much less contact with 
the physics of the problem than the dynamic approach 
adopted here. What work did exist on the dynamics, in- 
cluding our own, has errors or inconsistencies which we 
hope that we have rectified. What is required most of 
all are more experiments in order to guide the theoret- 
ical development, ruling out and supporting the various 
theoretical approaches. We hope that the current work 
stimulates the carrying out of such experiments in the 
future. 



Acknowledgments 

RGM wishes to thank the EPSRC (UK) for the award 
of a postgraduate studentship. 



[1] P. L. Luisi, The Emergence of Life (Cambridge Univer- 
sity Press, Cambridge, 2006), Chap. 10. 

[2] U. Seifert, Adv. Phys. 46, 13 (1997). 

[3] W. Helfrich, Z. Naturforsch. C 28, 693 (1973). 

[4] Z-C. Ou-Yang, L. Ji-Xing and X. Yu-Zhang. Geometric 
Methods in the Elastic Theory of Membranes in Liquid 
Crystal Phases. (World Scientific, Singapore, 1999). 

[5] D. Fanelli and A. J. McKane, Phys. Rev. E78, 051406 

(2008) . 

[6] S. R. de Groot and P. Mazur. N on- Equilibrium Thermo- 
dynamics (Dover, New York, 1984). 

[7] R. G. Morris and A. J. McKane. In preparation. 

[8] B. Bozic and S. Svetina, Eur. Biophys. J. 33, 565 (2004). 

[9] B. Bozic and S. Svetina, Eur. Phys. J. E24, 79 (2007). 
[10] B. Bozic and S. Svetina, Phys. Rev. E80, 013401 (2009). 
[11] D. Fanelli and A. J. McKane, Phys. Rev. E80, 013402 

(2009) . 

[12] Z-C. Ou-Yang and W. Helfrich, Phys. Rev. Lett. 59, 2486 
(1987); Z-C. Ou-Yang and W. Helfrich, Phys. Rev. A 39, 



5280 (1989). 

[13] C. J. Adkins. Equilibrium Thermodynamics. (Cambridge 
University Press, Cambridge, 1983). Third edition. 

[14] O. Kedem and A. Katchalsky, Biochim. Biophys. Acta 
27, 229 (1958). 

[15] O. Kedem and A. Katchalsky, Trans. Faraday Soc. 59, 
1918 (1963). 

[16] W. H. Beyer. CRC Standard Mathematical Tables (CRC 

Press, Baton Rouge, 1987). 
[17] http://www.wolfram.com 

[18] M-P. Nieh, J. Pencer, J. Katsaras, and X. Qi, Langmuir 

22, 11028 (2006). 
[19] L. Beney, E. Linares, E. Ferret, and P. Gervais, Eur. 

Biophys. J. 27, 567 (1998). 
[20] J. Pencer, G. F. White, and F. R. Hallett, Biophys. J. 

81, 2716 (2001). 
[21] L. Miao, U. Seifert, M. Wortis, and H-G. Dobereiner, 

Phys. Rev. E49, 5389 (1994). 



